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I.  INTRODUCTION 


The  volume  of  fibrous  composites  used  as  high  performance  structural  materials  is  increasing 
very  rapidly  because  of  the  possibilities  to  adapt  the  mechanical  properties  of  the  composites  in  the 
elastic  domain  to  a  particular  loading  situation.  This  tailoring  of  the  extcnsional  and  bending 
stillnesses  as  well  as  of  the  coupling  between  bending  and  extension  is  not  possible  for  the  fatigue 
life  or  the  static  strength.  The  failure  modes  of  composite  materials  are  many  and  interdependent. 
Consequently,  the  understanding  of  the  phenomenon  has  to  be  related  to  the  manner  in  which  li¬ 
bers  and  matrix  which  have  different  mechanical  properties,  are  put  together  to  form  a 
unidirectional  composite. 

Many  attempts  have  been  made  in  the  past  to  describe  the  properties  of  a  composite  lamina 
based  on  the  properties  of  its  constituents.  In  all  of  these  studies  hypotheses  have  been  made  about 
the  stress  transfer  between  the  fibers  and  the  matrix.  Tliis  load  transfer  is  the  key  to  the  microme¬ 
chanical  approach  of  composite  materials.  The  properties  of  both  the  matrix  and  the  fibers  can 
be  determined  experimentally  without  difficulty,  but  if  the  load  transfer  is  not  understood,  the  state 
of  stress  in  the  different  constituents  cannot  be  evaluated,  and  consequently,  the  properties  of  the 
composite  cannot  be  calculated.  Until  now,  the  mathematical  representations  of  the  stress  transfer 
between  the  fibers  and  the  matrix  have  given  results  that  do  not  compare  well  with  experimental 
measurements. 


I.  INTRODUCTION 


I 


Much  work  has  been  done,  especially  by  chemists,  to  define  and  analyze  the  zone,  known  as 
the  fiber  matrix  interface,  where  the  stress  transfer  occurs.  It  has  been  proved  with  chemical  analysis 
and  photomicrographs  [lj  that  a  layer  of  finite  thickness,  exists  between  the  fiber  and  the  matrix, 
and  that  it  is  chemically  different  from  both  of  them.  This  layer,  described  in  references  2  and  3, 
is  referred  to  as  the  interphase,  and  poses  a  new  problem  because  its  thickness  is  only  on  the  order 
of  a  micron.  As  a  result,  the  mechanical  properties  of  this  region  cannot  yet  be  measured  exper¬ 
imentally. 

It  is  now  commonly  accepted  that  progress  in  micromechanics  is  not  possible  without  both 
an  understanding  of  the  physics  of  the  interphase  and  its  mathematical  representation.  Conse¬ 
quently,  interdisciplinary  studies  (4,5|  have  determined  that  the  mechanical  simulation  of  the  stress 
transfer  between  the  matrix  and  the  fibers  must  include  parameters  that  result  from  a  better 
understanding  of  the  chemistry  of  adhesion.  According  to  this  idea,  a  mechanical  model  of  stress 
transfer  that  includes  an  interphase  layer  has  been  developed  [6],  and  this  study  improves  the  pre¬ 
vious  analysis  by  considering  the  possibilities  of  damage  in  the  interphase. 

The  fiber  matrix  stress  transfer  is  not  the  only  problem  in  the  determination  of  the  properties 
of  unidirectional  composites.  The  interaction  between  fibers  has  been  observed  to  be  of  significant 
influence  when  discontinuities  exist  in  the  filaments  [7j.  Tiber  discontinuities  can  be  attributed  to 
the  use  of  whiskers,  or,  in  continuous  fiber  composites,  to  flaws  in  the  filaments  and  stresses  in¬ 
duced  by  the  fabrication  and  curing  processes.  During  the  service  life  other  fiber  breaks  appear  due 
to  the  loading  of  the  structure.  Using  the  results  for  the  stress  concentrations  in  a  three  dimensional 
array  of  fibers  reported  in  reference  7,  this  present  study  examines  the  influence  of  the  fiber-fiber 
interaction  on  the  stress  transfer  between  the  matrix  and  the  fibers. 
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II.  LITERATURE  REVIEW 


Our  study  is  oriented  towards  developing  a  mathematical  model  of  the  stress  transfer  between 
the  matrix  and  the  fibers.  However,  as  mentioned  in  a  later  section,  different  models  correspond 
to  the  different  types  of  bonding  between  the  fibers  and  the  matrix  that  have  been  observed  in 
composite  materials.  Consequently,  for  the  micromechanical  study  of  fibrous  composites  to  be 
significant,  the  physics  of  the  bonding  process  needs  to  be  understood. 


A.  Chemistry  of  Fiber-Matrix  Bonds 


Along  with  the  understanding  of  the  fiber-matrix  bonds,  another,  by  no  means  insignificant, 
purpose  of  studying  the  chemistry  of  adhesion  is  the  improvement  of  the  bonds. 

Many  studies  have  already  provided  results  about  glass/epoxy  systems  [4-6],  the  composites 
used  in  the  experimental  part  of  this  research.  These  systems,  because  of  the  very  high  strength  of 
the  fiber  matrix  bond  in  dry  environment,  are  often  used  for  data  comparison  in  the  literature. 
However,  this  report  is  not  limited  to  glass/epoxy  composites.  In  other  cases  where  no  chemical 
link  is  established  between  the  fibers  and  the  matrix,  the  interphase  is  only  a  thin  layer  of  the  matrix 
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that  has  reacted  differently  from  the  bulk  material  because  of  the  presence  of  the  liber  This  layer 
is  not  mechanically  attached  to  the  liber  and  therefore  constitutes  a  weak  boundary  layer.  The  at¬ 
tachment  of  the  interphase  to  the  fiber  can  be  realized  in  most  of  the  cases  by  applying  different 
surface  treatments  to  the  fibers,  or  by  adding  coupling  agents  to  the  matrix.  Hence,  chemical 
analysis  techniques  have  replaced  the  empirical  scries  of  tests  that  have  been  used  in  the  past  to 
investigate  fiber-matrix  bonds.  In  reference  8,  Baun  reviews  no  less  than  fifty  four  microscopic, 
spectroscopic,  thermodynamic,  and  kinetic  techniques  of  investigation  of  the  interphasc. 

The  chemistry  of  the  matrix  is  usually  available,  but  the  surface  of  the  fibers  has  to  be  ana¬ 
lyzed.  Glass  surfaces  have  been  studied  using  different  techniques  listed  by  Sabat  in  reference  |6|. 
Furthermore,  in  order  to  obtain  a  better  understanding  of  the  functionality  of  the  surface,  a  deri¬ 
vation  process  adapted  from  Everhart's  and  Reilly's  technique  |9]  can  be  used.  This  procedure 
consists  of  using  reagents  that  contain  elemental  tags.  These  reagents  are  chosen  to  react  with 
specific  functional  groups,  and  by  scanning  for  an  elemental  tag  using  XPS,  the  presence  of  the 
corresponding  functional  group  can  be  identified  on  the  surface. 

It  must  be  said  that  all  the  techniques  of  analysis  of  surface  chemistry  cannot  be  used  once 
the  fibers  and  the  matrix  are  put  together.  Consequently,  the  nature  of  the  interface  inside  actual 
composites  remains  a  guess  based  on  the  analysis  of  the  elements  before  assembly. 


B.  Mechanics  of  the  Fiber- Matrix  Adhesion 


The  mechanics  of  matrix-fiber  load  transfer  has  been  studied  since  the  early  1950  s  1 1 0],  but 
it  was  not  until  196.1  that  the  stress  concentration  due  to  discontinuities  in  the  fibers  was  considered 
|  Ilf  Along  with  the  analytical  research,  many  experimental  studies  have  been  conducted  on  the 
subject  adapting  photoelasticity  and  other  test  measurements. 
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1.  Analytical  Approach 


Originally  the  problem  of  the  matrix-fiber  stress  transfer  was  approached  by  applying  fiber 
theories  1 10-  14J.  Later,  four  different  types  of  analysis  arose; 

-  Elasticity  theory 

-  Finite  element  analysis 

-  Fracture  mechanics 

-  Shear-lag  analysis 

These  are  referenced  in  detail  by  Sabat  (6|  who  also  classifies  the  papers  on  the  subject  according 
to  the  assumptions  made  about  the  bond  and  the  nature  of  the  interface  between  the  fibers  and  the 
matrix. 

The  treatment  of  the  stress  concentration  due  to  discontinuities  in  the  fibers  has  been  studied 
separately  [7],  and  has  rarely  been  included  in  mathematical  models  of  the  matrix  fiber  stress 
transfer.  The  only  cases  when  the  influence  of  the  neighboring  fibers  has  been  taken  into  account 
are  statistical  treatments  of  the  unidirectional  composite  strength  115,16]. 


2.  Experimental  Approach 


The  measurement  of  the  maximum  stress  transfer  between  the  matrix  and  the  fibers  can  be 
approached  in  two  ways.  Mechanical  tests  such  as  short  beam  shear,  four  point  shear,  or  flexural 
strength  can  be  run  on  composite  specimens  to  determine  indirectly  the  Interfacial"  shear  strength. 
But  the  results  of  these  tests  are  very  controversial  because  of  the  complex  states  of  stress  they  in¬ 
duce  at  the  interface  between  the  matrix  and  the  fibers. 

A  more  direct  approach  is  to  measure  the  maximum  stress  transfer  by  observing  the  behavior 
of  a  single  fiber.  Four  different  single  fiber  tests  have  been  developed  thus  far: 
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-  The  pull  out  test.  The  experiment  is  represented  in  figure  la),  and  the  following  test  procedure 
is  used:  a  single  fiber  is  embedded  in  a  thin  disc  of  the  polymer  matrix  in  such  a  manner  that  it 
represents  the  axis  of  the  sample.  Then  the  fiber  is  pulled  out  of  the  polymer,  and  the  necessary 
force  is  recorded.  The  thickness  of  the  polymer  disc  is  increased  until  the  liber  cannot  be  pulled 
out  and  breaks.  The  pull  out  force  is  plotted  versus  the  thickness  of  the  disc,  which  is  also  the 
length  of  the  bond.  The  slope  of  this  curve  is  supposed  to  be  the  "interfacial  shear  strength"  |17|. 

The  high  scatter  in  the  data  led  the  researchers  to  question  the  assumption  of  pure  uniform 
shear  stress  at  the  interface  ( 1 8, 19].  It  appears  that  the  polymer  meniscus  at  the  reentrant  corner 
generates  stress  concentration,  and  that  a  tensile  stress  exists  that  pulls  the  matrix  away  from  the 
fiber.  Consequently,  variations  in  the  shape  of  the  polymer  block  have  been  investigated  |20).  l  he 
pull  out  droplet  type  specimen  is  represented  in  Figure  lb).  Because  the  analysis  of  this  test  is  still 
not  precise  about  what  is  measured,  the  results  of  this  experimental  procedure  must  be  considered 
as  being  qualitative  information. 


-  The  microdebond  test  (F  igure  2).  1'his  is  the  most  recent  test  method,  and  probably  the  most 
representative  of  the  problem  we  are  concerned  with  because  it  studies  a  fiber  embedded  in  an  actual 
composite  material.  A  cross  section  of  a  fibrous  composite  is  polished,  and  placed  under  a  micro¬ 
scope.  A  spherical  indenter  is  then  placed  over  the  end  of  a  fiber  and  loaded  until  the  detachment 
between  the  fiber  and  the  matrix  occurs.  The  load  recorded  at  the  time  of  the  detachment  senes 
as  input  in  a  micromechanical  model  that  gives  the  "interfacial  strength."  Rcgretably,  the  under¬ 
standing  of  the  mechanism  of  detachment  is  not  complete,  and  there  is  no  evidence  that  the  re¬ 
corded  load  is  a  function  of  the  debonding  process  only  |21,22). 


-  Transverse  tensile  test  |23|.  A  single  fiber  is  embedded  in  a  compression  specimen  of  the 
polymeric  matrix.  F  'ndcr  load,  the  Poisson  s  effect  creates  a  tensile  force  normal  to  the  fiber  in  the 
center  of  the  specimen.  The  debonding  of  the  fiber  from  the  matrix  is  detected  using  reflected  light, 
and  the  corresponding  value  of  the  load  is  used  to  calculate  an  "intcrfaciai  tensile  strength." 
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figure  2 


Schematic  of  the  microdebond  test 
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-  Critical  fiber  length  test.  Because  it  is  easy  to  run,  this  test  method  has  beer  largely  used  and 
studied,  l'he  initial  fiber  length  test  was  originally  described  by  Kelley  |24|.  A  latter  section  of  this 
report  contains  a  detailed  description  of  the  test,  and  proposes  some  modifications  in  order  to  tailor 
it  to  the  study  of  fibrous  composite  materials. 
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A.  INTRODUCTION 


The  mechanical  properties  of  fiber  reinforced  composites  are  significantly  influenced  by  the 
stress  transfer  between  the  matrix  and  the  fibers.  Several  approaches  arc  used  to  quantitatively 
approximate  this  load  transfer.  Rosen,  followed  by  many  others,  considers  the  liber-matrix  bond 
as  being  perfect  [13].  This  means  that  the  fiber  and  the  matrix  constitute  a  single  unit  with  two 
zones  having  different  mechanical  properties.  Furthermore,  the  'dual"  material  behaves  elastically, 
even  in  highly  stressed  regions  like  the  fiber  tips.  This  model  has  been  greatly  improved  by  I*.  J. 
Sabat  who  adds  an  interfacial  layer  between  the  fiber  and  the  matrix.  But  still,  all  the  materials  aie 
considered  as  elastic,  and  the  possibilities  of  breakage  or  yielding  of  the  matrix,  or  of  slippage  of  the 
fiber  with  respect  to  the  matrix  are  not  considered. 

The  purpose  of  the  first  model  developed  in  this  chapter  is  to  quantify  the  possibilities  of 
damage  in  an  "elasto-damageable"  matrix  surrounding  a  fiber  in  the  case  of  a  strong  chemical  bond, 
i.e.,  in  the  case  of  continuity  at  the  atomic  level  between  the  fiber  and  the  interphasc,  and  between 
the  interphase  and  the  matrix.  This  study  is  conducted  for  a  fiber  embedded  in  a  fibrous  composite 
material,  and  considers  all  the  problems  of  stress  concentration  due  to  broken  fibers.  In  the  second 
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model,  the  fiber  matrix  adhesion  is  assumed  to  be  only  mechanical  with  no  atomic  link  between  the 
fiber  and  the  matrix.  Both  models  are  then  used  to  define  the  ineffective  length.  The  influence  of 
the  ineffective  length  on  the  mechanical  properties  of  short  fiber  composites  has  been  addressed  in 
the  past  by  Reifsnidcr  et  al.  [47] .  and  consequently  is  not  reiterated  in  this  report. 

i 

j  B.  SHEAR  LAG  MODEL  WITH  A.V  INTERFACIAL 

!  LA  YER  (SINGLE  FIBER) 


In  the  last  decade  much  progress  has  been  made  in  the  field  of  adhesion  science  One  of  the 
most  important  advancements  is  obviously  the  observation  in  certain  cases  of  an  interfacial  layer 
that  exists  between  the  adhesive  and  the  adherend.  As  a  result,  P.  J.  Sabat  developed  a  shear-lag 
model  |6|,  derived  from  Rosen  s  theory  |l.f|,  to  calculate  tell  stresses  in  the  fiber,  the  neighboring 
matnx  and  the  interfacial  layer,  for  a  single  fiber  embedded  in  a  composite  material.  I  lis  hypotheses 
are: 

•  the  fibers  carry  only  tensile  load 

'  •  the  matrix  carries  only  shear  stress 

•  the  average  material  carries  only  tensile  load 

•  the  interphase  carries  part  of  the  tensile  load  of  the  fiber  and  part  of  the  shear  stress  oi 
the  matrix 

fr,(.v)  -  ,1  •  (7;( r)  (!) 

Tjix)  =  li  .  t„a  x)  (0 


•  shear  transfer  is  limited  to  matrix  interphase 
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•  the  bonds  between  fiber-interphase,  interphase-matrix  and  matrix-average  material  . im¬ 
perfect. 

The  hypotheses  of  the  shear-lag  analysis  have  already  been  used  by  many  researchers  i>' 
ported  by  Holister  and  I'hoinas  [28],  and  are  reasonable  enough  in  the  case  of  high  modulus  libers 
embedded  in  a  ductile  matrix.  However,  the  last  assumption  concerning  the  perfect  bond,  even  near 
the  fiber  end,  has  been  criticized  by  Amirbayat  and  Iiearle  [29)  who  consider  the  possibility  of 
fiber-matrix  debonding  and  slippage  at  the  fiber  end  because  of  the  high  stress  concentration  in  this 
region. 

Sabat  also  neglected  the  overstress  in  the  matrix  due  to  discontinuities  in  the  neighboring  li¬ 
bers.  Because  he  was  more  interested  in  comparing  resin  properties,  filament  surface  characteristics 
and  surface  treatments,  this  approximation  was  not  significant  for  his  results.  This  simplification, 
however,  certainly  cannot  be  used  to  characterize  the  stress  field  in  a  fiber  embedded  in  a  composite 
material. 

Using  Sabat's  equation  for  the  equilibrium  of  the  unit  cell  represented  with  our  coordinate 
system  in  Fig.  4,  we  obtain  the  following  formula  for  the  stresses  in  the  fiber  and  in  the  neighboring 
matrix:  (the  resolution  for  the  change  of  axis  is  given  in  Appendix  1) 

a  Ef 

Oj(x)  =  -p— ( 1  -  exp[  -  i;x])  (3) 

and 

-  Rfr\  a  l.f 

Tm(x)  - - — -  exp[  -  vx]  (4) 

where: 

2  2  -  Rl)  +  2piEfRr  _ 

( Rm  -  Ri)IGm  +  fl{Rt  -  Rfj/Gi 
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This  approach  of  the  fiber-matrix  adhesion  could  be  the  storting  point  for  a  new  generation 
of  models  which  take  into  account  the  existence  of  an  interfacial  layer  whose  chemical  and  me¬ 
chanical  properties  are  different  from  those  of  both  adhesive  and  adherend.  The  influence  of  the 
interphase  upon  the  strength  of  the  material  itself  will  be  given  in  a  later  part  of  this  report. 
However,  for  equations  3  and  4  to  be  representative  of  the  conditions  in  a  fiber  reinforced  plastic, 
the  overstresses  due  to  discontinuities  in  the  neighboring  fibers  have  to  be  taken  into  account. 
Furthermore,  this  model  has  shown  that  perfect  fiber-matrix  bonding  is  probably  an  illusion  near 
the  fiber  end.  Consequently,  modifications  have  to  be  made  in  order  to  consider  this  phenomenon. 


C.  STRESS  CONCENTRA  TIO N  SEAR 
DISCONTINUITIES  OF  FIBERS 


The  problem  of  stress  transfer  between  a  broken  fiber  and  the  neighboring  ones  has  already 
been  extensively  studied,  and  different  types  of  2-D  models  have  been  d,  .eloped:  linear  elastic- 
analysis  [30,31),  elastic  plastic  analysis  [32|  and  bilinear  stress  strain  curves  (33)  have  been  included 
in  finite  element  programs;  analytical  solutions  have  also  been  found  using  the  shear  lag  theory 
(7,34,35).  The  most  advanced  model,  to  our  knowledge,  is  that  of  Hedgepeth  and  Van  Dickc  ( 7) 
and  presents  the  solution  of  the  3-D  problem  for  fibers  arranged  in  an  hexagonal  array,  Figure  5, 
using  a  shear  lag  analysis. 

Fhe  figures  in  Table  1  have  been  calculated  using  tills  last  reference,  and  the  following  hy¬ 
potheses: 

•  Zero  to  six  of  the  outer  fibers  can  break. 

•  The  fiber  breaks  are  infinitely  small, 
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•  The  breaks  of  the  outer  fiber  induce  stress  concentration  at  the  interface  between  the 


central  fiber  and  the  matrix. 

•  In  the  plane  perpendicular  to  the  fibers  and  passing  tlirough  the  break,  the  stress  con¬ 
centration  decreases  when  the  distance  from  the  break  increases. 

•  The  fiber  volume  fraction  is  equal  to  .44. 

•  The  value  of  the  stress  concentration  induced  in  the  matrix  by  a  number  of  coplanar 
breaks  of  the  outer  fibers  is  constant  around  the  central  fiber,  and  it  is  equal  to  an  average 
value  calculated  at  the  center  of  the  fiber. 

We  notice  that  the  first  value  corresponds  with  the  result  of  the  finite  element  analysis  con¬ 
ducted  by  Barker  [30], 

A  point  neglected  in  almost  all  the  previous  treatments  is  that  the  overstress  is  not  only  ap¬ 
plied  to  a  plane  but  to  a  volume  and  this  in  a  non-equal  fashion.  The  distribution  of  overstress 
has  been  studied  in  both  two  dimensions  [27,36,37,38]  and  three  dimensions  [7],  We  consider  in 
this  study  that  the  stress  concentration  is  maximum  in  the  plane  of  the  break,  and  that  it  decreases 
along  the  axis  of  the  central  fiber  proportionally  to  the  distance  from  tliis  plane.  Wc  assume  that 
the  effects  of  coplanar  breaks  of  the  outer  fibers  is  not  sensitive  at  the  interface  between  the  central 
fiber  and  the  matrix  if  the  distance  from  the  plane  of  the  breaks  is  larger  than  eight  times  the  radius 
of  the  fibers.  This  value  of  the  length  of  the  overstressed  zone  has  been  chosen  according  to  pre¬ 
vious  finite  element  studies  of  the  problem  [30,  33) .  Consequently,  the  mathematical  function  re¬ 
presenting  the  overstress  in  the  matrix  along  the  interface  between  the  central  fiber  and  the  matrix 
is  defined  by: 


•6) 


and  is  shown  in  figure  6. 
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Figure  6.  Variation  of  shear  stress  concentration  near  an  infinitely  small  discontinuity 
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D.  PERFECT  BOND  (MULTIPLE  FIBERS) 


1.  Effects  of  the  Stress  Concentration  on  the  Shear  Lag  Model 


The  stress  concentration  has  so  far  been  defined  in  the  matrix  only;  hence,  the  shear  stress  in 
the  matrix  rm  is  formulated  as  follows: 

For  0  <  x  <  Si R, 


and  for  x  >  8 Rf 


—  a  Rflij  Y] 

T^(.v)  = - -  - —  cxp[  -  IJ.vJ 


Because  the  equilibrium  of  the  fiber  still  has  to  be  satisfied,  the  equation  7A  of  Appendix  1  is  still 
valid  and  is  recalled  here 


-  2/? 

~dF  ~  rJ~  Tm 


By  integrating  equation  9  we  obtain  for  0  x  < 


cxp(  --  >5 


v)x(r'  +  xk/1  ~r‘])JrfY^ 


where  K  is  a  constant  of  integration. 
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Performing  the  integration: 


zIIl 

t- 


1  ~  Cl  tl.X  +  1  ,  r  ,, 

ci  +  ^ —  )  txp[  -  „,]  +  A  | 


1  he  boundary  condition  <r/o)  =  0  inserted  in  Eq.  10  leads  to 


- c/  + 


1  - 


SR, 


HI) 


Eq.  1 1  into  Eq.  10 


oj(x)  = 


]  —  c • 

c,(  exp[  -  rjx]  -  1)  +  '(b*  +  1]  exp[  -  *ix~\  -  1) 

V>KfY] 


(12) 


Figures  7  and  8  show  respectively  the  distribution  of  the  shear  stress  in  the  matrix  and  the  tensile 
stress  in  the  fiber. 

It  is  useful  to  remember  that  all  the  calculations  in  this  section  have  been  made  for  infinitisely 
strong  linear-elastic  materials  .  As  we  shall  see  later,  this  is  far  from  reality  and  the  breakage  of  the 
fibers  occurs  at  a  stress  level  lower  than  the  one  for  which  the  discontinuities  in  the  stress-strain 
curves  appear.  Consequently,  Figures  8  and  9  are  displayed  here  uniquely  to  demonstrate  the  in¬ 
fluence  of  the  stress  concentration  due  to  discontinuities  in  the  neighboring  fibers  and  are  of  no 
value  concerning  the  actual  stress  field  in  and  around  a  fiber  embedded  in  a  fibrous  composite 
material  subjected  to  tensile  load. 
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Figure  8.  Elasticity  model:  tensile  stress  in  the  fiber  (for  one  to  six  broken  neighboring  fibers) 
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2.  Ineffective  Length 


a.  Plastic  Behavior 


Rosin  i!  reference  -  definr*  the  inctTectiu*  length  is  1<  lliv.w 

\l  ji'ino  di-la!. ve  Iri'tn  ,:;i  internal  t'r‘..ik.  li'vr  stress  will  bt  .1  given  : 
lil'cr  stress  ~  f  >■  1  One  mav  define  (his  fraction  O  ol  the  average  si.. 

■  over  vdii.h  (lie-  sue--  i-  I.--  i!',ni  a  a  )  hi.iv  be  con-ldend  iia  !ii  ■.  :i,yj 
ki.Ctl.  •"  is  deli'ii. .!  : . 


(/l  --  (T,(f1).''<T,|C-o) 


1  he  mode!  with  an  intertaei.it  laser  developed  in  Appendix  1  results  i 


./>  —  I  c\p[  •  - 


this  leads  to: 


d  =  1  jt]  x  In [  1  /( 1  -  d>) J 


Inserting  l;q.  5  into  1  q.  15  we  obtain: 


<5  =  1- 


( Rm  —  Rt)l(im  +  /h/\,-  -  R,V(,, 


2Rt!W  ~  K)  *  2^>  Rt 


'n f  1(1  --  <i> 


Assuming  a  perfect  bond,  i.e  .  a  good  stress  ti  in-tc  f .  ■;  me 

assuming  the  elastic  behavior  of  the  materials.  a  li!  >:■  .-  ;  A-  ' 

stress,  ('onsequcntly,  the  possibility  of  hrnk.u’e  in  'h  •  avtliv  to. i . - : ru ' 1 1 
1  he  eMcvt'  ot  f'v.  nrpn.o  ,l.t  •  ; 


tonmila  ol  the  -tress  m  if."  li!--.:  n  .r 


i|i.  >.\  ..t  v  1 1<  M  W  1  '-i-. 
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I  i  s  [(»;<’  *-  it  e\p( 


<  (  1  -  exp[  -  »;  vj »  -  i<  - 


Hits  equation  i<  solved  numetwalk  to  Ini. I  <\  mil  !  icuv  g  «!n'  .  ■  t 


ccntr.it:>'::  eh  :.:.:cn  cd  b\  <  .  n  •!,  *  :<  ;■  :.i 

j 

-< 


h.  Hduirior  in  the  i'aniavci!  State 


In  tl'C  pt.-'. fits  section,  the  behavim  ut  the  liber.  matrix.  rt:>i  nr  ::  ’  e  -uuu  J  :>•  he 

elastic.  even  though  the  poh nier  matrix  is  known  to  \ :e;J  under  cd'.'...  h  :  -he. a  i« -..Jjttu  In 
this  >e.t:.  :t  the  hb:  : -s 'ti'l  eonsidneri  ei.e -;e.  hut  t-n!1!  t:..  !'■  J  :  .  •  i-mmed 

perlcdly  elastic -plastic.  A  schemas  •  of  the  behavior  of  the--,  vuitiriA  •  ■;  ite  tie  !'h 

Consequently  the  iu'm  of  tin*  shear  strength  of  the  m.i'nx  to  ::  w  t  :r,t. : phase 

n.  part  deteuninc'  ’he  locus  of  failure  Similarly.  t.ic  t.uie  ot  tile  .a  yield::  p. ri  the  matrix 
to  the  shear  \ lelding  point  o!  the  intcrpha.se  Jctem  lines  v. hv’i.er  the  tn  itrix  or  the  mt. 'phase 
umlc:  goes  plasiie  deformation. 

i  lie  following  mathematical  mode!  ..an  ..a  •■>. : . ’ - .  .  i  s  !iv'cr-t.uti;\  s-. 51  ym  w  ah  chemical  links 
between  the  constituents.  But  in  order  to  use  tin.  model  in  an  appiopn  t:e  tnannet .  the  composite 
s\  stem  must  he  cliaracteri/.ed  according  to  whether  the  debonding  occurs  because  of  the  failure  of 
the  n*-.-r phase  nr  of  the  failure  ol  the  matrix  1  lie  locus  ol  the  failure  can  he  determined  b\  using 
surf  i.  ■:  1:1  i!> •  is  techniques  on  the  surface  of  the  fiber  as  described  in  reference  h  |  hen,  if  the  failure 
",  .nrs  in  the  interphase,  our  model  can  he  used  to  cal*  white  cither  the  ineffective  length  or  i!  e  bond 
' t r .. 1 1 ig 1 1 : .  ...•  cr.ling  to  which  of  these  is  used  as  input  data  On  the  other  band,  it  the  failure  occurs 
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the  case  of  the  glass  epoxy  composite  used  in  the  experimental  pan  of  tics  re'eauh  I  lie  me¬ 
chanical  properties  of  the  constituents  are  listed  in  Table  2. 

As  shown  in  1  mure  7.  the  matrix  evidences  a  large  damaged  /one.  and  men  a  region  where 
rupti'R  occur'  near  the  fiber  tip.  However,  whether  'lie  behavior  o'  the  m  itux  i>  cln  he  or  :vt.  Site 
same  load  has  ;<•  be  translated,  t onsv q  tenth.  the  atu  under  the'  two  cur  e  in  I  mure  1  1  have  to 
be  equal 

1  o  define  the  function  representing  the  shear  stress  in  the  in  ,t..  damaged  state,  the 

length  Xr.  X  1  and  Xd  nut'1,  he  defined.  1  irst  of  all.  the  leorth  Xt  a1'  n  wh  ,c  rep  fur :  occurs, 
can  be  defined  by. 

•  :  •  !  s ...  (IS) 

where  S„  is  the  shear  strength  of  the  matrix. 

I  he  mathematical  function  represented  in  1  igutc  7,  and  used  in  this  section,  to  characterize 
the  overstress  in  the  matrix  near  the  fiber  end  shows  a  singular  point  at  a  distance  from  the  fiber 
tip  equal  to  eight  times  the  radius  of  the  fiber.  As  a  result,  we  now  have  to  solve  the  equations 
where  the  o  vet  stress  function  appears  in  two  steps  in  order  to  consider  the  two  different  domains 
of  definition  of  the  function. 

Thus  tw  o  formula  can  be  used  for  r'(.r)  w  hether  the  value  of  the  shear  stress  at  the  end  of  the 
overstressed  zone  r„,(8 R.)  is  smaller  or  larger  than  S„.  the  shear  strength  of  the  matrix. 

Consequently,  for  t„(S R.)  >  (the  ease  represented  in  figure  8),  fq.  ?  is  insetted  m  1  q 
1 S 

~  5  Kf  f'-,  ’?  (  v>  \ 

■  "~:/7  J; - txp[  -  »7  A/]  X  (  C,  +  ~~  [  I  -  c,  ] )  -  S„.  <1‘)| 

1  he  *.  i -!•  t  w  •  'o  -''Ive  tins  cqu  item  is  numerical!'.  ficur  12  ie| "••••••  'Us  ]  '  1  ot  ■  V 

'ticss  ,  ■  ••x •,  nt :  in  -u  cm  die  icneth  Xr 

ft',,  the  it-.-.t  ...me.,  b'r  V  .>  r _ ■ 2\ . , .  1  q  s  is  inserted  in  1  q  IS 

m.  vs  m  i  t  tv  vi  \i  \  -is  :a 


Table  2.  Input  data  for  the  elasto-damaged  model 


& 


a  average  material 
m  =  matrix  material 


i  =  interphase 
f  =  fiber 


Ef  =  71  800  MPa 
Rf  =  62.5E-06  m 
Xf  =  410.0  MPa 

•  Gi  =  15  000  MPa 

Ri  =  63.SE-06  m 

Gm  =  765  MPa 
Rtn  =  80.7E-06  m 
Tym  =  21.37  MPa 

c  Sm  =  32.0  MPa 

Ea  =  2  700  MPa 
Ra  =  IE-03  m 

0-1. 

•  a  =  15.44  MPa 
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RfEf* 


2/5  i:a 


exp[  -  tf  Xr]  =  Sn 


(20) 


* 


or 


Xr- -  l/v  ln[  -2P  Ea  SJa  RfEfii]  (21) 

In  both  cases,  it  is  assumed  that  no  load  is  transferred  between  the  matrix  and  the  fiber  from  the 
end  of  the  fiber  to  the  point  of  abscissa  Xr. 

The  length  Xd  of  the  zone  where  damage  occurs  must  now  be  defined.  The  solution  can  be 
found  through  the  formula 


rXI 

t m(x)dx  =  r y(Xd  -  Xr)  ( 22) 

Jo 

which  equates  the  surfaces  under  the  two  curves  in  Figure  11,  and  where  ry  is  the  yielding  shear 
stress  of  the  matrix. 

In  order  to  solve  Eq.  22  we  need  to  define  XI  the  point  for  which  the  shear  stress  takes  the 
value  of  the  yielding  point  of  the  matrix.  The  equation  for  the  shear  stress  in  terms  of  XI  is  given 
by: 


Tm(*^0  ~  Ty 

The  shear  stress  is  represented  by  either  Eq.  7  or  Eq.  8  whether  the  value  of  the  shear  stress  at  the 
end  of  the  overstressed  zone  is  smaller  or  bigger  than  ry  the  shear  yielding  point  of  the  matrix. 
•  Consequently,  as  it  has  been  done  before  to  solve  Eq.  18,  Eq.  23  is  solved  in  two  steps.  Eq.  22  is 

first  solved  for  the  case  when  t m(ftR/)  >  ry,  and  Eq.  7  is  inserted  into  Eq.  23: 


a  i:fRfy} 
'2/?  La  ' 


exp[  —  rfXf] 


C  + 


XI 

HRf 


<24 1 
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*V»  y,j’  *  >  '  j»  *",>  v 


VI 


This  equation,  as  Eq.  19,  is  to  be  solved  numerically,  and  Figure  13  gives  the  solution  as  a  function 
of  the  stress  concentration  factor  c,.  In  this  case,  when  r„(.x)  for  x  between  zero  and  XI  is  entirely 
described  by  Eq.  7,  Eq.  22  becomes 


—  a  ErRiY] 


exp[  -  r,xV 


1  —  <y)^  dx  =  t  y(\d  —  Xr) 


—  a  Ef  Rfy\  rxi  1  -  r ,  rxi 

— 2p~f~ -  exp[  -  yjx]dx  +  j  xexp[  -  rjx]dx  =  ry(.\  d  -  X r) 


i-o  r*' 


Xd=  'ia  lRf  0<exp[  -  VXI]  -  1)  +  ±-^L((r1Xl+  1)  exp[  -  r,Xt]  -  1)1  +  Xr  (27) 

ZP  ba  ry  [_  8/</ 


where  XI  and  Xr  are  respectively  given  in  Eq.  24  and  Eq.  19.  The  influence  of  the  stress  concen¬ 
tration  on  the  length  of  the  damaged  zone  is  represented  in  Figure  14. 


For  rm(8 Rf)  <  -ry  Eq.  8  is  inserted  into  Eq.  23: 


XI  =  —  !/>;  x  ln[  -2/?  Ea  Tyjo  EfRj\\] 


In  this  case,  t„(x)  for  x  between  zero  and  XI  is  described  partly  by  Eq.  7,  and  partly  by  Eq.  8.  As 
a  result,  Eq.  22  becomes 
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Figure  13.  Length  XI 
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-  -  - .. 


/V 


—  a  EfRjti 


cxp[  -  n*]/  c(  +  (1  -  c,) \  dx 


—  a  EfRrti 

+  —2PE  . .  CXPC  "  =  Xy{Xd  ~  Xr) 


and  by  integrating  Eq.  29  we  get: 


o  E  R  I 

xd=  cK  expc  _  8/?^]  ~l)  + 1^7  ((8 R /*»  +  0  «pC  -  %Rfri  ~  0 


+  exp[  -  rfXt]  —  expf 


After  the  point  of  abscissa  Xd,  the  matrix  stays  in  its  elastic  domain  and  the  shear  stress  function 
t£(jc)  needs  to  be  defined  in  order  to  calculate  the  tensile  stress  af(x)  in  the  fiber.  In  order  to  ac¬ 
complish  this,  we  note  that,  for  0  <  X  <  Xr 


4(*)  =  « 


a}(x)  =  0 


For  Xr  <  x  <  Xd 


4(*)  =  fy 


Inserting  Eq.  33  into  Eq.  9  and  integrating  we  get: 


°f(x)  =  — Ty(x  -  Xr) 
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1 W  1  '.F  '  VF  '  ^ '  W  W  rj  r:rnrjra  w\i  w-* 


~M  w  X  **  1* 


For  Xd  <  X,  the  shear  stress  r^fx)  in  the  matrix  can  be  obtained  by  a  simple  translation  of  the 
function  represented  cither  by  Fq.  7  or  by  Eq.  8.  The  choice  between  Eq.  7  or  Eq.  8  depends  on 
whether  the  length  of  the  damaged  zone  in  the  matrix  is  smaller  or  larger  than  eight  times  the  radius 
of  the  fiber,  the  length  of  the  overstressed  zone.  The  amplitude  of  the  translation  is  equal  to  the 
difference  of  the  length  of  the  damaged  zone  to  the  legnth  XI. 

For  Xd  <  x  <  8 Rf  +  Xd  -  XI,  and  recalling  Eq.  7: 

Tm(*)  =  — TTr~~  exp£  ~  7(*  -Xd+  XI)]  x 
t-a 

Integrating  Eq.  9. 


.  x  -  Xd  +  XI  , , 

Cf+ — s*— (1_Ci) 


(35) 


odf(X)  =  -  2  mf 


and  inserting  Eq.  35  into  Eq.  36  we  obtain: 


fXd  r* 

dx  +  Tm(x)dx 
JX>  JXd 


(36) 


ad/{x)  = 


-ip  Ty  a  EfY)  r  -  c, 

(Xd  -  Xr)  +  —d—  cxp[  -  nix  -  Xd  +  AT)] 

exp[  -  n(x  ~  Xd  +  A/)]] 


ry 

I  -  y  nix  -  Xd  +  XI) 


8  R 


V 


Jxu 


and  developing  gives: 


-2  ^ 

a  fix)  = 


R 


V  ^  Ey* 

y  iXd-Xr) - -L 


V 


Ci  + 


1  —  Cj  n(x  -  Xd  +  A  Y)  +  1 
~M 2 


x  exp[  -  nix  -  Xd  +  XI)]  -  (  q  +  ~  '  )  expf  -  *;A7 ] 


(57) 


For  Xd  <  8 Rf  +  Xd  •  XI  <  X,  the  shear  stress  r£(x)  can  be  represented  by  the  translation  of  l  q 
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(38) 


d  ,  V 

*«(*)  = 


—  a  RfEjt\ 

~JX~ 


exp[  -  t\(x  —  Xd  +  AT)] 


Then,  the  integration  of  Eq.  9  gives: 


*/<*>-- f- 


/, 


JM  r»R,+Xd-Xl  . 

Tydx+\  rd 

Xr  JXd 


dm{x)dx  +  [ 


ZRf+Xd-Xl 


(39) 


Inserting  both  Eq.  35  and  Eq.  38  respectively  in  the  second  and  third  integrals  of  Eq.  3°  we  get: 


a  Ef 


od/(x)  =  -f-ry(xct-Xr)-  — 


1  -  ct  ffXl+l 


-\c,+ 


8  Rf  1 


C‘  +  \^L  x  exp[  -  SRftjl 


exp[  —  tfXf]  +  exp[  -  ij(x  -  Xd  +  A 7)]  -  exp[  -  8fy>j] 


,(40) 


For  8 Rf  +  Xd  ~  Xl <  Xd<  X,  rdJx)  is  represented  by  the  translation  of  Eq.  8  as  in  Eq.  38,  and  in¬ 
tegrating  Eq.  9,  we  obtain: 


dx  + 


x—Xd+XI 


r 

Jxt 


■'mix)dx 


(41) 


and  inserting  Eq.  38  into  Eq.  41  we  get: 


»/(*)  =  T y  (xd  -  xr)  -  [exp  [  -  r)(x  -Xd+  A 7)]  -  exp[  -  AT]]  (42) 


All  the  logical  cases  have  been  reviewed  and  are  summarized  in  a  flow  chart  represented  in 
Figure  15.  The  overstressed  line  follows  the  path  used  to  generate  the  curves  in  Figures  13  and  14 
after  the  data  reported  in  Table  2.  Figures  16  and  17  show  the  influence  of  the  stress  concentration 
on  the  shear  stress  in  the  matrix  and  the  tensile  stress  in  the  fiber,  respectively. 
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Figure  15.  Flow  chart  of  the  computer  program  to  compute  the  ineffective  length 
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Figure  17.  Elasto-damaged  model:  tensile  stress  in  the  fiber  (for  one  to  si*  broken  neighboring  fibers) 
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3 
3 
m, 


The  tensile  stress  in  the  fiber  is  now  totally  defined,  and  Rosen  s  formula  lor  the  fiber  eili- 
ciency  can  be  used  again  to  fmd  the  ineffective  length  S.  In  figure  IS,  the  applied  stress  is  chosen 
in  order  to  obtain  in  the  center  part  of  the  fiber  a  stress  intensity  equal  to  the  fiber  strength  I  his 
maximum  value  is  a  material  characteristic,  and  thus,  the  ineffective  length  still  appears  as  a  func¬ 
tion  of  the  stress  concentration  only. 


3.  Results  and  Discusssion 


1  he  present  analytical  model  improves  the  accuracy  of  the  mathematical  representation  of  the 
matrix-fiber  stress  transfer  by  considering  the  possibility  of  damage  in  the  matrix  or  in  the  niter- 
phase  and  the  stress  concentration  near  fiber  breaks.  As  a  result,  the  equations  which  describe  the 
stress  fields  in  the  matrix  and  in  the  fiber  take  complex  forms,  thus  the  particular  influence  of  the 
different  parameters  is  difficult  to  evaluate. 

The  influence  of  the  stress  concentration  around  liber  breaks  is  illustrated  in  figure  IS.  I  his 
curve  exhibits  a  minimum  value  which  is  characteristic  of  the  liber-inatnx  system.  I  he  location 
of  this  minimum  on  the  stress  concentration  axis  is  related  to  the  ratio  of  the  fiber  strength  to  the 
yielding  point  of  the  matrix.  One  must  also  notice  in  f  igure  IS  that  when  the  stress  concentration 
increases,  the  ineffective  length  does  not  drop  but  remains  in  the  same  order  of  magnitude.  I  his 
is  not  the  case  when  the  calculations  are  performed  only  in  the  elastic  domain.  If  this  tendency 
were  verified  experimentally,  single  fiber  samples  could  be  used,  at  least  for  a  first  approximation, 
to  evaluate  the  efficiency  of  the  matiix-fiber  stress  transfer  in  composites. 

Experimental  measurements  of  the  ineffective  length  are  available  in  the  literature  |4A4n|, 
and,  as  described  later  in  this  report,  we  have  conducted  several  series  of  experiments  of  our  own 
and  taken  many  measurements,  for  the  comparison  to  be  possible  be'ween  experimental  meas¬ 
urements  and  analytical  values  obtained  through  the  elasto-damaged  model,  a  short  parametric 
study  must  be  conducted  to  examine  the  influence  of  the  main  material  properties  figures  10  and 
20  show  the  influence  of  the  radius  of  the  fiber  and  of  the  shear  modulus  of  the  matrix  respective!) 
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on  the  ineffective  length.  In  order  to  refer  Ln  f  igure  20  to  the  different  tspes  ot  resin--  uwd  in  the 
previous  studies,  we  also  varied  the  values  of  the  shear  yielding  point  and  of  the  shear  strength  ot 
the  matrix.  The  mechanical  responses  under  shear  loading  of  the  different  matrices  used  in  I  lgtire 
20  are  reported  in  Table  3. 


E.  LIMITED  FRICTION 


1.  Stress  Field  in  the  Fiber 


For  many  fiber-matrix  systems,  the  model  previously  developed,  assuming  the  perfect  bond¬ 
ing  of  the  fibers  to  the  matrix,  gives  much  lower  results  than  the  experimental  measurements  of  the 
ineffective  length  with  polarized  light.  For  thermoplastic  matrix  composite  materials  by  example, 
the  polymer  does  not  undergo  further  polymerization  during  the  fabrication  of  the  composite. 
Consequently,  it  is  less  likely  to  react  with  the  fibers  than  a  thermosetting  resin  The  same  ahscencc 
of  chemical  bonding  can  also  occur  in  glass/cpoxy  systems  when  the  fibers  undergo  certain  surface 
treaments.  Also  in  the  case  where  there  is  chemical  bonding  between  the  matrix  and  the  fibers,  the 
interface  may  constitute  a  weak  boundary  layer  which  breaks  before  the  matrix  reaches  its  damaged 
state. 

In  these  two  cases,  the  possibility  of  fiber  slippage  in  the  matrix  reported  in  reference  4  must 
be  considered.  Consequently,  in  this  section,  a  mathematical  model  using  the  basic  theory  of 
elasticity  is  developed  for  the  shear  stress  at  the  interface  and  for  the  tensile  stress  in  the  liber 

For  fibrous  composite  materials  in  which  the  fibers  arc  not  bonded  to  the  matrix,  and  w  hich 
;ire  loaded  in  tension,  the  shear  stress  in  the  matrix  is  generated  by  frictional  forces  due  to  the  latent! 
contraction  of  the  matrix.  T  his  can  be  represented  by: 


J5 
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Ef  7.18x10 
Rf  6.25xl0“8 


I.5XI010 
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Table  3.  Mechanical  response  of  the  matrices  used  to  generate  Figure  20 


Gm 

(MPa) 

t  yielding 
(MPa) 

Sm 

(MPa) 

550 

15.36 

25.3 

650 

18.15 

29.9 

750 

21.37 

34.5 

850 

23.74 

39.1 

950 

26.53 

43.7 
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Ti  =  n  Pj  ( l ) 

where  r ,  is  the  shear  stress  in  the  matrix  at  the  interface,  n  is  the  coefficient  of  friction  between  the 
fiber  and  the  matrix,  and  P,  is  the  lateral  pressure  at  the  interface. 

The  solution  for  equation  1  is  found  tlirough  the  equilibrium  of  the  unit  cell  represented  in  figure 
21.  The  matrix  and  the  fiber  are  both  considered  isotropic  linear  elastic  materials.  Consequently, 
the  principle  of  superposition  can  be  applied  as  follows: 

Equilibrium  of  the  unit  cell  in  Figure  21: 

n  Rf  af+n(Rm~  Rf)°m  =  I'  ^ 


Equilibrium  of  the  fiber  (Figure  21): 


2n  Rj-Xj  dx  +  n  Rf 


d  Of 

-d-dx  =  0 
dx 


After  simplification  of  Eq.  3: 


Rf  d  Of 

r,s=~~dT 


(3) 


(-D 


The  pressure  P ,  can  be  obtained  by  considering  the  compatibility  of  the  radial  deformation  at  the 
interface: 


i’/F  +  UfPi  -  UmF+  (5) 

where  Ufr  is  the  radial  deformation  of  the  fiber  due  to  the  applied  load.  ( 'fP]  is  the  radial  deformation 
of  the  fiber  due  to  the  applied  load,  UfP  is  the  radial  deformation  of  the  liber  due  to  interfacial 
pressure,  VmT  is  the  radial  deformation  of  the  matrix  due  to  the  applied  load,  and  l.’mr  is  the  radial 
deformation  of  the  matrix  due  to  interfacial  pressure. 
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Determination  of  vfF 

The  fiber  is  assumed  to  undergo  only  tensile  load  as  represented  in  Figure  23a.  We  write 
Hooke's  law  for  the  fiber: 


<7=c/e* 


and  by  definition  of  Poisson's  ratio 


Eq.  7  is  introduced  into  Eq.  6: 


The  integration  of  Eq.  8  gives: 


-  VfOf 

V^—U-R, 

V 


Determination  of  UfPj 

The  fiber  is  now  loaded  in  lateral  compression  only  as  represented  in  Figure  23b.  Equilibrium 
of  the  fiber  can  be  written  as  follows: 


7r  —  —  Pj  and  erg  =  —  Pj 


We  now  write  Hooke's  law  for  the  fiber: 


Cf,  =  ~  v/<rr)  =  £r 


And  the  strain  displacement  relationsliip  reads: 


III.  ANALYTICAL  ANALYSIS 


liWNTi, 


i.n.r 


c 


A  r  =  r 


ze 


Introducing  Eq.  10  and  1 1  into  Eq.  12  we  obtain: 


P,Rf 
uf?r— 


Determination  of  Umf 

Eq.  2  is  used  to  define  <r„  the  tensile  stress  in  the  matrix  represented  in  Eigure  23c 


F  —  it  Rj  oj 


<Rl  -  R}) 


Erom  Lame  theory  |39] 


Er  -  tg  ~  tH 


Consequently  the  strain  displacement  can  read: 


Ar  =  ~vmar-  vm  ° h) 


Boundary  conditions: 


aH  =  am  <  °e(Rf)  =  0  ;  ar(Rf)  =  0 


Including  these  three  boundary  conditions  into  Eq.  16  gives: 
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Determination  of  Um  r 

As  before,  I amc  theory  is  used  to  solve  the  problem  shown  in  Eigure  23d  of  the  matrix 
undergoing  the  interfacial  pressure  only. 
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Wc  write  Hooke's  law  lor  the  matrix: 


Ar  =  -p— (<*//  -  vm  <rf) 


l  IS) 


The  equilibrium  of  the  matrix  gives: 


°H  =  P, 


(RJr)2  +  1 
(*«/*/  -  1 


(19) 


and 


<V 


(RJr)2  -  1 

(RJR/  -  1  _ 


(20) 


For  r  =  R/,  Eq.  20  becomes: 


(21) 


and  Eq.  19  becomes: 


(*«/*/)  +  1 
(PJKj)  -  l 


(22) 


Eq.  21  and  22  into  Eq.  18  give: 


P,Rf\  (/Uty2  +  i 

L  (KJKf)2  '  > 


(23) 


All  the  displacements  having  been  defined  by  Eq.  9,  13,  17,  and  23,  arc  now  included  into  the 
compatibility  equation,  Eq.  5. 
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or,  introducing  the  constants  5,  S3,  a,  b,  and  c,  Eq.  24  becomes: 


P,  =  SxF-S2of 


where 


S,  =  1  jn(R2m  -  Rj) 


a1  +  1 
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and 


a 2  +  1 
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and 


*/ 


Introducing  Eq.  25  into  Eq.  1 


t(  =  m[5,  f-Sjay] 


Eiq.  28  into  Eq.  4: 


Rj-  drsj 

T“5T 


=  m[5,  F- S2  ay-] 


and  developing  Eq.  29,  we  obtain 


c/<r 


dx 


7  ,  _2n 


>2ff/  = 


2m 


Because  the  end  of  the  fiber  is  free  of  stress,  the  solution  for  Eq.  30  is 
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Figure  24  shows  the  repartition  of  the  tensile  stress  along  the  fiber  for  a  unit  force  I7  and  for  different 
values  of  the  friction  coefficient. 
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3.  Results  and  Discussion 


It  must  be  noted  that  the  formulation  of  this  model  is  very  basic,  and  is  presented  in  this  re¬ 
port  only  in  order  to  complement  the  elasto-damaged  model.  A  more  elaborate  model  of  limited 
friction  or  a  model  such  as  Adams  (41]  could  be  applied  to  the  length  where  the  rupture  of  either 
the  matrix  or  the  interface  occurs. 

Figure  25  shows  the  influence  of  the  coefficient  of  friction  p  on  the  ineffective  length.  Fq. 
34  illustrates  that  the  value  of  6  depends  only  on  the  material  properties.  W  hen  slippage  occurs, 
the  ineffective  length,  as  just  defined,  has  the  same  significance  as  in  the  case  of  perfect  bond,  as  tar 
as  he  stress  transfer  between  the  matrix  and  the  fibers  is  concerned.  However,  we  have  not  found 
a  good  solution  to  verify  experimentally  the  results  given  by  Fq.  34.  1  he  measurement,  under 

polarized  light,  of  the  deformation  zone  in  the  matrix  becomes  difficult  because  of  fiber  slippage. 

This  model  of  limited  faction  does  not  take  mto  account  any  interfacial  layer  because  the 
interphase  has  been  defined  as  a  weak  boundary  layer,  and  thus  does  not  constitute  a  structural 
bond  as  in  the  "perfect  bond"  model.  The  contribution  of  the  interphase  to  the  stillness  of  the 
composite  is  also  neglected  because  of  the  very  small  volume  of  this  layer  compared  to  the  other 
constituents. 

This  model  is  also  affected  by  the  mathematical  limitation  which  is  due  to  the  logarithmic 
fonn  of  the  solution  in  Tq.  35.  The  term  in  the  logarithm  cannot  be  negative,  and  as  a  result  the 
mechanical  and  geometneal  properties  of  the  materials  cannot  be  varied  freely  For  example,  the 
modulus  of  the  average  material  in  the  fiber-matrix  system  used  to  draw  Figure  25  cannot  be  less 
than  2.3E10  Pa.  Consequently,  the  single  fiber  test  where  the  average  material  is  the  resin  itself, 
cannot  be  investigated  with  this  model.  Another  problem  attributed  to  the  mathematics  of  the 
solution  is  the  impossibility  of  applying  this  model  to  a  case  without  slippage.  The  friction  coeffi¬ 
cient  is  in  the  denominator  of  the  fraction,  consequently  the  infinite  value  corresponding  to  a 
'perfect  bond"  yields  a  value  for  the  ineffective  length  equal  to  zero.  This  has  been  proved  false 
tlirough  both  experimental  anu  analy  tical  studies  of  the  matrix-fiber  stress  transfer. 
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IV.  EXPERIMENTAL  PROCEDURES  AND 


OBSERVATIONS 


A.  DESCRIPTION  OF  THE  MA  TERIALS 


Fibers 


Three  types  of  glass  fiber  have  been  used  in  this  study: 

•  Coming  optical  fibers  (outer  diamter  125  ^m;  core  diameter  85  nm) 

•  Coming  optical  fibers:  The  original  outer  diameter  of  the  fibers  is  125  (im. 
The  filaments  arc  etched  in  a  bath  of  hydrofluoric  acid  to  obtain  a  diameter  of 
80  /im.  The  diameter  of  the  fibers  decreases  at  the  rate  of  3  /im/min.  Thus  the 
fibers  stay  in  the  bath  of  acid  during  15  minutes,  and  then  are  rinsed  with  dis¬ 
tilled  water. 

•  Vetrotex  Saint-Gobain  E-glass  fibers  (diameter  24  /im) 

We  used  an  epoxy  system  composed  of  Fpon  828  resin  produced  by  Shell  and  cured 
with  either  35,  45,  or  50  weight  parts  of  jeffamine  D230  produced  by  Texaco 
Chemicals.  The  curing  cycle  of  the  specimens  is  described  in  Figure  26. 
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Matrix 


interphase  The  acrylate  coating  is  removed  from  the  optical  libers  with  a  razor  blade.  All  the 
fibers  are  first  cleaned  with  acetone  then  thoroughly  rinsed  with  distilled  water. 
Some  of  the  fibers  are  etched  in  hydrofluoric  acid. 


B.  PROPOSAL  FOR  A  FIBER-MATRIX  ADHESION 


ADVANCED  TEST  METHOD 


1.  Background 


The  fiber  matrix  appears  to  be  a  determining  factor  of  the  fibrous  composite  materials' 
strength.  As  long  as  the  stress  transfer  between  the  matrix  and  the  fibers  is  not  quantified,  the 
micromechanical  approach  of  the  mechanical  properties  of  composites  will  remain  an  inapplicable 
theory. 

A  method  largely  used  in  the  past  to  quantify  the  stress  transfer  between  the  matrix  and  the 
fibers  consists  of  embedding  a  single  fiber  in  a  thermoplastic  or  thermoset  matrix  and  then  sub¬ 
mitting  the  matrix  to  a  tensile  load  in  the  direction  of  the  liber  |24-29|.  If  the  fiber  matrix  bond  is 
strong  enough  to  transfer  a  load  able  to  break  the  filament,  a  minimum  lcgntli  of  the  liber  pieces 
will  be  reached,  Figure  27.  In  the  particular  cases  of  a  single  fiber  embedded  in  a  transparent  matrix, 
the  critical  fiber  length  can  be  measured  with  an  optical  microscope.  If  the  matrix  is  opaque,  it  can 
be  dissolved  and  filtered  to  separate  the  small  pieces  of  fiber 

The  critical  fiber  length,  lc,  is  characteristic  of  the  maximum  load  that  can  be  translered  from 
the  matrix  to  the  fiber  by  shear  Kelley  derived  a  simple  expression  for  the  shear  stress  at  the 
interface: 
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where  oc  is  the  tensile  stress  in  the  fiber,  and  d  is  the  diameter  of  the  fiber.  As  explained  by  fra/.er 
[25,26|  this  value  corresponds  to  the  shear  strength  of  the  weakest  link  which  can  be  either  the 
matrix  near  the  fiber  or  the  interfacial  layer  between  the  fiber  and  the  matrix.  Consequently,  it  is 
more  accurate  to  refer  to  tc  as  the  stress  transfer  coefficient". 

1  his  test  can  be  modified  to  study  the  effects  of  the  fiber-liber  interaction  on  the  fiber  matrix 
stress  transfer.  To  do  so,  several  fibers  have  to  be  embedded  parallel  and  close  together  in  the  same 
specimen.  Bascom  and  Jensen  [27|  have  been  the  pioneers  in  this  domain,  using  specimens  that 
contain  two  to  five  fiber  aligned  in  the  same  plane.  In  keeping  with  the  study  of  fibrous  composites, 
a  sample  containing  a  three-dimensional  array  of  libers  is  developed  in  this  study. 


2.  Description  of  the  Sample 


As  shown  in  f  igure  4,  the  dogbone  shape  of  the  single  fiber  specimen  is  kept,  but  seven  fibers 
are  now  embedded  in  the  matrix  in  a  hexagonal  configuration. 

In  order  to  take  measurements  of  the  critical  fiber  length  on  the  central  fiber,  it  has  to  be 
possible  to  differentiate  it  from  the  others.  Consequently,  we  have  used  optical  fibers  without 
plastic  coating,  embedded  in  a  transparent  epoxy  resin.  Thus,  besides  carrying  mechanical  load,  the 
glass  fibers  provide  an  optical  path.  A  laser  light  is  conducted  only  through  the  central  fiber,  and 
the  specimen  is  loaded  in  tension.  As  described  before  the  fiber  breaks  in  several  parts,  and  at  each 
of  the  breaks  a  part  of  the  light  intensity  is  transmitted  through  the  break  to  the  other  piece  of  the 
liber,  and  the  rest  of  the  light  is  dispersed  in  the  surrounding  matrix  signaling  the  location  of  the 
break  It  the  proportion  of  light  transmitted  to  the  other  part  of  the  fiber  is  large,  it  is  possible  to 
sec  several  breaks  We  have  been  able  during  our  experiments  to  observe  a  maximum  of  twenty 
three  breaks  I  he  light  transmission  through  the  breaks  has  been  considerably  improved  by  the 
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usage  of  multi-mode  fibers.  These  fibers  have  a  bigger  diameter  core  and  consequently  a  larger 
acceptance  cone.  The  other  problem  encountered  in  the  study  of  epoxy  specimens  is  the  very  low 
strain  to  failure  of  the  matrix  itself  This  inconvenience  has  been  avoided  with  the  help  of  the 
Virginia  Tech  Polymers  l  aboratories  which  developed  a  system  which  undergoes  seven  percent 
elongation  before  failure. 

The  production  of  these  specimens  is  tedious,  and  performed  manualv.  Tins  procedure  does 
not  allow  the  use  of  small  diameter  fibers,  and  as  a  result,  does  not  allow  the  study  of  the  types  of 
fibers  normally  used  in  composite  materials.  Consequently,  a  device  to  organize  automatically  the 
fibers  in  the  mold  has  to  be  developed  to  remove  this  limitation  before  this  test  method  can  be  of 
any  utility. 

C.  SPECIMENS  PREPARATION 

This  study  requires  four  different  types  of  specimens: 

-  NFAT  Rl.SIN  DOGBONF.  SPFCIMF.NS:  They  are  machined  from  a  plate  obtained  by  casting 
bulk  resin  in  the  following  manner,  degassed  matrix  is  poured  into  an  open  air  mold  that  is  kept 
perfectly  horizontal  in  order  to  obtain  a  thickness  as  constant  as  possible  with  the  lowest  residual 
stress  due  to  the  casting.  After  machining,  the  specimens  arc  grinded  and  dry  polished.  The  ge¬ 
ometry  of  the  specimen  is  shown  in  Figure  28. 

-  IQSPIFSCC  SPFCIMFNS:  These  are  machined  from  the  same  plate  as  the  dogbone  specimens. 
1  he  geometry  of  the  specimens  is  shown  in  Figure  29. 

-  SINGFF  F1BFR  SPFCIMFNS:  Single  fiber  specimens  are  obtained  by  following  the  process 
related  by  P.  .1  Sabat  in  reference  6.  The  fibers  arc  first  mounted  in  the  mold  as  indicated  in  Figure 
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30.  and  then  degassed  matrix  is  poured  into  the  mold  at  room  temperature.  The  TRV  3032  silicone 
rubber  produced  by  Dow  Coming  has  been  chosen  to  make  the  mold  because  of  the  particular 
properties  of  this  material:  accurate  reproduction  of  standard  dogbone  specimens,  nonadhesion  to 
epoxy  resin,  high  temperature  resistance,  and  easy  usage.  The  final  geometry  of  the  samples  is 
shown  in  Figure  30. 

-  SEVEN  FIBER  SPECIMENS:  A  bundle  of  seven  fibers  is  first  prepared  out  of  the  mold. 

Step  1:  Seven  fibers  are  placed  into  a  plastic  tube,  its  inner  diameter  equal  to  three  times  the  ex¬ 
ternal  diameter  of  the  fibers.  The  fibers  arrange  themselves  in  an  hexagonal  array  according  to  the 
dimension  of  the  tube.  A  six  inch  strand  of  each  fiber  is  left  out  of  the  tube.  The  bundle  is  kept 
in  position  in  the  tube  by  a  droplet  of  glue. 

Step  2:  The  part  of  the  fibers  outside  the  tube  is  stripped  of  its  coating  with  a  razor  blade,  dipped 
into  acetone  to  remove  the  remaining  particles  of  coating,  rinsed  with  distilled  water,  and  in  certain 
cases,  etched  with  NaOlI. 

Step  3:  1  he  fibers  are  placed  on  a  silicone  rubber  plate  and  maintained  in  position  with  two  pieces 
of  tape  as  shown  in  Figure  31.  The  free  ends  of  the  fibers  are  then  bonded  together  with  a  droplet 
of  glue. 

Step  4:  The  center  parts  of  the  fibers  are  brought  closer  to  each  other  by  tightening  two  pieces  of 
very  thin  thread  around  them  as  shown  in  Figure  31. 

Step  5:  The  bundle  is  mounted  m  the  same  mold  used  to  prepare  the  single  fiber  specimens,  and 
the  degassed  matnx  is  poured  into  the  mold 
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Figure  31.  Steps  3  and  4  of  the  7  fiber  specimen  preparation 
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D.  MECHANICAL  TESTING 


1.  Tensile  Testing  of  the  Fibers 


In  order  to  obtain  the  tensile  strength  of  the  fibers  as  a  function  of  its  length,  a  tensile  test  is 
performed  according  to  AS  I'M  3379-75  on  an  Instron  macine.  model  1  122  s  n  4470.  I  he  test  is 
run  for  three  different  gage  lengths  with  a  constant  cross  head  speed  of  0,2  in  mn.  1  he  results  arc 
presented  in  Table  4. 


2.  Shear  Testing  of  the  Matrix 


We  demonstrate  in  the  analytical  part  of  this  report  that  the  shear  modulus,  the  shear  yielding 
point,  and  the  shear  strength  of  the  matrix  have  a  crucial  influence  on  the  approximation  of  the 
ineffective  length,  and  consequently  on  the  evaluation  of  the  strength  of  the  composite  itself. 

A'  a  result,  the  shear  test  methods  commonly  used  have  been  reviewed  to  select  the  most 
accurate  with  regard  to  the  measurement  of  the  three  mechanical  properties  we  are  interested  in. 
f  rom  a  theoretical  point  of  view,  the  most  uniform  shear  stress  is  obtained  by  applying  a  torsional 
load  to  a  thin  wall  tube,  but  this  test  cannot  determine  the  shear  strength  because  failure  occurs  due 
to  the  buckling  id’  the  specimen  Previous  studies  |42]  have  demonstrated  that  the  torsion  of  a  solid 
rod  is  not  the  solution  cither,  because  of  the  nonlinear  distribution  of  the  shear  stress  along  the 
radius  of  the  rod  in  the  case  ot  plastic  deformation.  The  plate  twist  method  and  the  split  ring  test 
weie  abo  abandoned  because  they  measure  only  the  shear  modulus.  Similarly,  the  short  beam 
shear  test  is  not  suitable  became  it  measures  only  the  shear  strength  Other  test  methods  such  as 
the  tail  shear  test  and  the  cross  sandwich  beam  test  are  not  considered  because  stress  concentrations 
m  the  tested  pail  of  the  specimens  make  their  results  questionable 
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Table  4.  Tensile  properties  of  the  fibers 


Radius  in  pm 

12 

40 

f>2.5 

Gage  length  in  mm 

200 

120 

50 

Hf  in  MPa 

65000 

72100 

71S(K) 

Strength  in  MPa 

802 

701 

650 

Weibuil  parameters 


3.59 

2628.7 


3.09 
1 127.3 


f  inally,  the  simplest  test  satisfying  all  the  requirements  is  known  as  the  Iosipeseu  shear  test 
from  the  name  of  Nieolae  Iosipeseu,  who  lirst  proposed  this  test  in  the  beginning  of  the  at 

the  I'niversity  of  Bueharest.  The  test  was  primarily  designed  to  study  metals,  but  it  has  been  largely 
used  to  test  materials  such  as  adhesives  and  adhesive  bonds  in  1  ockhccd  laboratories,  and  com¬ 
posite  materials  at  NASA  l.angley  and  VPl&Sld  (42.43.44).  The  different  methods  used  to  inves¬ 
tigate  this  test  are  listed  in  Reference  3  and  the  conclusions  can  be  summarized  as  follows: 

-  Lite  Iosipeseu  specimen  produces  a  zone  of  uniform  shear  stres  large  enough  to  perform  meas¬ 
urements  (shaded  area  in  figure  2l>). 

-  The  notches  have  to  be  rounded  to  avoid  stress  concentrations  that  would  cause  the  failure  of  the 
specimen  before  the  actual  shear  strength  ot  the  material  is  reached. 

I  he  mechanical  analysis  of  the  test  is  given  in  the  original  paper  [45j  by  Nieolae  Iosipeseu  and  is 
confirmed  in  reference  42. 

Test  description: 

The  tests  arc  run  in  the  NASA  Program  laboratories  at  YPI&SU  using  the  fixture  designed  by 
D.  b.  Adams  and  shown  in  F  figure  32  with  a  specimen  mounted.  The  tests  arc  performed  with  a 
cross  head  speed  of  0.02  inches  per  minute. 

Because  of  the  error  introduced  by  the  compressive  deformation  at  the  landing  points,  meas¬ 
uring  the  relative  displacement  of  the  two  halves  of  the  fixture  is  not  rut  accurate  way  to  obtain  the 
shear  strain  at  a  point  located  on  the  line  between  the  notches'  tips.  Consequently,  wc  use  rec¬ 
tangular  rosettes  that  we  bonded  between  the  notches  to  measure  the  shear  strain.  The  Micro 
Measurements  rosettes  are  positioned  with  the  45  degree  gage  at  the  center  of  the  specimen  as 
shown  m  figure  33,  following  Walrath’s  and  Adams’  procedure,  'flic  smallest  strain  gages  have 
been  chosen  in  order  to  fit  into  the  zone  of  uniform  deformation.  According  to  Iosipeseu  himself 
and  also  to  later  photoelastic  studies  |42|  this  zone  has  a  width  ol  2.5  mm.  The  data  acquisition  is 
performed  through  the  MAI  PAC  software  developed  by  .1.  Ilidde  at  VPl&Sl,',  and  the  results  are 
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I'igure  32.  Iom'|h:nlu  lest  fixlmc  (sjiriifmn  inruiitt  il > 
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post  processed  with  Supercalc  4  to  correct  the  errors  due  to  transverse  sensitivity  and  gage  misa¬ 
lignment.  The  gage  misalignment  is  measured  under  optical  microscope  before  testing  the  sample. 

Six  tests  have  been  run  and  the  results  arc  reported  in  Table  5.  1  he  curve  in  Figure  34  has 
been  recorded  during  the  loading  of  the  specimen  number  one.  I  his  test  has  been  conducted  until 
specimen  failure  occurs  without  being  interrupted  by  the  buckling  of  the  sample.  This  has  not  been 
the  case  for  tests  numbers  two,  four,  and  six. 


3.  Observation  with  Polarized  Light  of  the  Single  Fiber  Sample 


The  epoxy  system  used  in  this  study  is  a  birefringent  material.  Consquently,  it  is  possible, 
using  polarized  light,  to  study  the  state  of  stress  of  the  matrix  near  a  liber  break.  The  accuracy  of 
tlus  experimental  technique  that  has  been  sued  with  success  in  the  past  [6,46),  has  been  improved 
in  this  research  by  the  usage  of  high  quality  optics  and  tensile  devices. 

The  specimen  is  first  loaded  in  tension  with  a  micro  tensile  rig,  the  MINIMA!',  built  by 
Polymer  laboratories  and  represented  in  Figure  35.  I'liis  apparatus,  small  enoguht  to  mount  on 
the  stage  of  a  microscope  and  controlled  by  computer,  provides  the  stress-strain  curve  of  the  tests. 
Once  loaded  to  its  maximum  stress,  the  specimen  is  either  kept  loaded  or  unloaded.  The  stressed 
zone  in  the  matrix  is  observed  through  a  polar  microscope  built  by  Spectra- l  ech.  Figure  36.  For 
each  of  the  fiber-matrix  systems,  fifty  breaks  are  investigated.  The  quality  of  the  optics  allows  us 
to  observe  different  zones  along  the  fiber,  as  shown  in  Plate  37: 


-  / one  1:  This  black  region  near  the  tip  of  the  fiber  is  probably  due  to  the  breakage  of  either  the 
matrix  or  the  interface.  The  size  of  this  zone  remains  the  same  when  the  specimen  is  unloaded. 


-  /one  2:  1  his  bright  write  area  is  the  stressed  region.  Its  length  diminishes  when  the  specimens 

arc  unloaded  Consequently,  we  tend  to  think  that  when  the  load  is  applied,  both  plastically  and 
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Table  5.  losipcscu  shear  test  results 


Specimen  number 

1 

2 

3 

4 

5 

6 

Number  of  part 
of  D230 

35 

35 

45 

45 

50 

50 

Gm  (MPa) 

950 

900 

765 

748 

720 

732 

t  Yielding  (MPa) 

.02%  oiTset 

26.53 

21.37 

18.20 

Shear  strength 
(Mpa) 


43.7 


34.5 


(MPa  ) 


rigure  36.  Spectra-tcih  micrtmnpc  (ruiniinat  mmr.U  Ij 
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elastically  deformed  zones  are  evidenced,  and  that  in  this  case,  the  length  of  this  /.one  is  the  inef¬ 


V 


V 


fective  length,  When  the  load  is  released,  otuy  the  plastically  deformed  /.one  appears. 

-  Zone  3:  This  zone  along  the  center  part  of  the  fiber  has  a  uniform  grey  color,  and  represents  the 
unstressed  area  of  the  matrix. 


4.  Critical  Fiber  Length  Test 


Encouraged  by  the  correlation  between  our  theoretical  results  and  P.  J.  Sabat's  experimental 
measurements,  we  use  this  test  again,  even  if  it  has  not  been  entirely  satisfetory  in  the  past  [6.27] 
The  critical  fiber  length  test  allows  an  easy  qualitative  comparison  between  different  fiber-matrix 
systems,  but  the  quantitative  results  of  this  method  are  still  impossible  to  relate  to  any  of  the 
matrix-fiber  stress  transfer  theories. 

However,  in  the  present  study,  we  use  this  test  to  investigate  qualitatively  the  difference  be¬ 
tween  the  single  and  the  seven  fiber  specimens.  The  single  fiber  samples  are  loaded  in  tension  with 
the  Minimat,  but  because  of  the  maximum  loading  force  of  two  hundred  newtons  of  this  device, 
the  seven  fiber  samples  are  loaded  on  an  Instron  machine  (model  1 122  S/N  4470).  A  force  of  270N 
has  to  be  applied  to  the  seven  fiber  specimens  to  reach  the  maximum  number  of  breaks  of  the 
central  fiber. 

A  laser  beam  generated  with  a  10  mW  laser  gun  is  shined  into  the  central  fiber  of  the  seven 
fiber  samples  to  take  measurements  of  the  broken  segments.  I  o  ensure  that  the  maximum  power 
enters  the  optical  fiber,  the  laser  beam  is  first  focalized  through  a  convergent  lens,  and  then,  the  end 
of  the  fiber  is  brought  to  the  focus  of  this  lens  with  a  stage  equipped  with  micromctric  displacement 
screws.  The  setup  is  represented  in  f  igure  3<S. 

The  measurements  are  taken  on  a  Zeiss  microscope  at  different  magnifications  according  to 
the  diameter  of  the  fibers.  1  or  the  24  ^m  diameter  fibers  we  use  a  160X  magnification  with  a  5 
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precision.  For  the  80  and  the  125  diameter  fibers  we  use  a  4UX  magnification  with  a  20 
pim  precision. 

5.  Results  and  Discussion 

a.  Critical  Fiber  Length  Test 

Single  Fiber  Samples 

The  results  of  the  series  of  tests  are  presented  in  fable  6.  The  epoxy  matrix  cured  with  35 
weight  parts  of  jefTamine  D230  appears  to  be  too  brittle  at  room  temperature,  so  the  first  break  of 
the  fiber  results  in  the  immediate  failure  of  the  specimen.  Consequently,  these  specimens  are  tested 
under  a  100W  lamp  to  elevate  their  temperature  to  about  60'C. 

As  expected,  the  diameter  of  the  fibers  influences  severely  the  critical  length.  On  the  other 
hand  the  influence  of  the  shear  modulus  of  the  matrix  is  totally  masked  by  the  broad  variability  of 
the  results.  The  dependency  of  the  critical  length  on  the  stmegth  and  geometry  of  the  fiber  results 
in  a  standard  deviation  of  the  measurements  equal  to  about  one  third  of  the  average  value,  Tliis 
is  large  enough  to  hide  the  expected  variation  of  twenty  percent  due  to  the  modulus  of  the  matrix. 

The  influence,  if  any,  of  etching  on  the  critical  fiber  length  appears  to  be  smaller  than  the 
standard  deviation  of  the  results.  Because  we  do  not  observe  any  major  change  whether  the  fiber 
is  etched  or  not,  we  would  tend  to  think  tha  for  this  fiber  matrix  system  the  interface  is  stronger 
than  the  matrix,  even  when  the  fiber  is  not  surface  treated. 

Seven  Fiber  Samples 

The  results  of  this  series  of  tests  are  reported  in  fable  7.  We  tried  to  perform  these  exper¬ 
iments  under  the  same  conditions  as  the  single  fiber  test,  but  less  tests  have  been  run  because  of  the 
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I  able  6.  C  ritical  fiber  length  (single  fiber  sample) 


Critical  Number  of  D240  parts  in  the  matrix 


mini =  1490 


mini  ~  1410 


mini  -  2100 


750 

07"! 

li-  2130 

mini  =  1980 
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complexity  involved  in  making  the  specimens.  However,  the  influence  of  the  fiber  diameter  is  still 
sensitive. 

Some  post  testing  observations  of  the  specimens  lead  us  to  think  that  the  failure  process  is 
different  in  the  seven  fiber  samples  than  in  the  single  fiber  one.  After  testing  the  seven  fiber  samples, 
we  ground  and  polished  the  four  sides  of  the  broken  dogbones  so  we  could  observe  the  fiber  breaks 
from  different  angles  under  a  microscope.  We  obtained  through  this  observation  two  significant 
results.  First,  a  break  in  the  central  fiber,  located  by  a  spot  of  light,  is  always  in  the  same  plane  as 
at  least  one  other  fiber  break.  This  would  tend  to  prove  that  the  matrix-fiber  stress  transfer  is 
'good,'  and  that  the  stress  concentration  near  a  fiber  break  is  sensitive  for  the  neighboring  fibers. 
The  second  observation  is  also  related  to  the  number  of  coplanar  breaks.  In  the  six  specimens  we 
observed,  we  have  not  found  more  than  three  coplanar  breaks.  This  is  probably  because  we  used 
fibers  of  a  large  diameter  that  can  carry  a  fraction  of  the  laod  which  is  not  negligible  with  respect 
to  the  strength  of  the  samples.  Consequently,  the  rupture  of  more  than  three  fibers  in  the  same 
plane  leads  to  the  immediate  failure  of  the  sample,  and  this  is  probably  lh<  explanation  for  the  liigli 
value  of  the  average  critical  length  recorded  in  tliis  series  of  tests. 

b.  Measurement  of  the  Deformed  Zone  in  the  Matrix  Under  Polarized  Light 

Table  8  presents  the  measurements  of  the  rupture  zone  recorded  on  the  single  fiber  samples. 
This  value  can  be  obtained  fairly  accurately  with  a  standard  deviation  of  about  one  fifth  of  the  liber 
diameter,  because  of  the  well  defined  bordeT  of  the  rupture  zone.  Each  of  the  figures  in  Table  8  is 
the  average  of  one  hundred  measurements  (fifty  breaks).  This  observation  has  a  high  degree  of 
reproductibility,  and  the  standard  deviation  is  smaller  than  the  precision  of  the  measure  itself. 
Consequently,  the  small  effect  of  the  variation  of  the  matrix'  mechanical  properties  can  be  recorded. 
As  expected,  the  length  of  the  rupture  zone  decreases  when  the  strength  of  the  matrix  increases. 

For  the  specimen  still  being  loaded  at  constant  deformation  in  the  Minimat,  the  value  re¬ 
corded  for  the  length  of  ihe  deformed  zone  depends  greatly  on  the  time  between  the  loading  and 
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Table  8.  Observation  under  polarized  light, 


Rupture 

Number  of  D230  parts  in 

the  matrix 

length 

in  ix  m 

35wp 

25C 

35wp 

60C 

45wp 

25C 

50wp 

25C 

F 

12  n  m 

32  nm 

35  nm 

40  fim 

40  fim 

b 

e 

r 

r 

a 

d 

40  ixm 

sample 

90  ^m 

i 

u 

s 

62.5  ^m 

n 

128  fim 

142  n m 

153  ;im 
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the  measurement.  A  measure  taken  24  hour  after  the  test,  while  the  sample  is  still  loaded,  is  the 
same  when  the  sample  is  unloaded.  Tliis  is  not  surprising,  because  of  the  viscoelastic  behavior  of 
the  matrix.  Consequently,  we  try  to  take  the  measurements  immediately  after  the  breaks  occur  in 
order  to  record  the  value  that  corresponds  to  both  the  elastic  and  the  plastic  deformation  of  the 
matrix,  and  as  a  result,  corresponds  to  the  ineffective  length  as  defmed  in  the  analytical  part  of  this 
report.  The  results  of  this  series  of  tests  are  reported  in  Table  9.  Because  of  the  shallow  contour 
of  the  zone,  the  error  on  this  measurement  is  about  one  time  the  diameter  of  the  fiber.  This 
repressents  a  precision  of  about  ten  percent  which  is  small  enough  not  to  mask  the  variation  of  the 
ineffective  length  w'ith  the  mechanical  properties  of  the  matrix. 

The  measurements  of  the  deformed  zone  24  hour  after  the  sample  has  been  unloaded  are  re¬ 
ported  in  Table  10.  These  values  are  characteristic  of  the  permanent  damage  that  occurs  in  the 
matrix  when  a  fiber  breaks.  The  length  of  tliis  damaged  zone  does  not  change  after  24  hours;  we 
spot  checked  the  specimens  two  weeks  after  the  tests,  and  the  same  values  were  recorded. 
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Table  9.  Observation  under  polarized  light.  Ineffective  length. 


Ineffective 

Number  of  D230  parts  in  the  matrix 

length 

* 

in  pm 

35wp 

25C 

35wp 

60C 

45wp 

25C 

50wp 

25C 

• 

12  ^m 

F 

i 

b 

182  fim 

200  nm 

246  f<m 

# 

e 

r 

40  n  m 
r  etched 

a 
d 

sample 

breaks 

626  nm 

> 

i 

u 

s  62.5 

tr 

782  )iva 

931  /im 

995  >im 

» 


e 


e 


e 
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V.  COMPARISON  BETWEEN  THE 


PREDICTIONS  OF  THE 
"ELASTO-DAMAGED"  MODEL  AND  THE 
OBSERVATIONS  UNDER  POLARIZED 
LIGHT 


The  approximate  agreement  between  the  analytical  predictions  of  the  elasto -damaged  model 
and  the  measurements  taken  under  polarized  light  is  demonstrated  in  Figures  39,  40,  and  41.  The 
theoretical  and  experimental  results  are  represented  by  solid  lines  and  point  markers,  respectively. 
However,  we  remain  critical  in  regard  to  our  model.  Because  a  mathematical  model  is  based  on 
assumptions,  we  discuss  in  this  section  the  discrepancies  as  well  as  the  resemblances  between  theory 
and  experiments. 

Figure  39  which  represents  the  length  Xr  of  the  rupture  zone  along  the  fiber,  is  probably  the 
easiest  of  the  three  graphs  to  interpret.  Theoretically,  only  the  shear  strength  of  the  matrix,  or  of 
the  bond,  influences  the  value  of  Xr.  The  assumption  that  the  matrix  and  the  interphase  are  per- 
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fectly  elasti-plastic  does  not  introduce  any  error  in  this  result.  Consequently,  the  experimental 
measurement  of  the  length  Xr  is  a  reliable  input  data  for  calculating  the  bond  strength.  However, 
Figure  39  shows  a  sensitive  difference  between  our  measurements  and  the  analytical  curve.  The 
fracture  of  the  bond  always  appears  to  be  longer  than  the  calculated  value.  A  possible  explanation 
of  this  phenomenon  is  that  the  plane  of  maximum  stress  concentration  passes  through  the  crack 
tip,  and  consequently,  moves  along  the  fiber  when  the  crack  propagates.  Furthermore,  the  stress 
concentration  at  the  instant  of  the  break  can  possibly  be  larger  than  the  value  calculated  according 
to  reference  7.  This  instantaneous  phenomenon  can  be  attributed  to  the  sudden  relaxation  of  an 
important  quantity  of  energy.  The  influence  of  this  "dynamic  overshoot"  has  been  addressed  by 
Hedgepeth  [34|,  and  the  results  are  on  the  order  of  1.2.  Such  an  increase  in  the  stress  concentration 
would  significantly  influence  the  length  Xr.  But  we  have  not  included  this  factor  in  our  math¬ 
ematical  model  because  the  lack  of  results  and  literature  on  the  subject  prevents  us  from  comparing 
values  obtained  through  different  methods.  This  instantaneous  phenomenon  is  also  very  vaguely 
defined  in  the  literature,  and  to  model  it  mathematically  would  require  a  study  of  its  own. 

Figures  40  and  41  represent  the  length  Xd  of  the  plastically  deformed  zone  in  the  matrix  and 
the  ineffective  length,  respectively.  These  two  graphs  show  a  very  encouraging  agreement  between 
our  theory  and  the  measurements.  The  slightly  higher  values  of  the  experimental  results  can  be 
attributed  to  the  length  Xr  which  appears,  as  explained  previously,  larger  than  the  value  expected. 
However,  such  a  close  agreement  (less  than  9.5  percent)  is  unexpected,  considering  all  the  as¬ 
sumptions  we  have  made  concerning  the  behavior  and  the  mechanical  properties  of  the  matrix  and 
of  the  interphase  in  particular.  As  a  result,  a  more  thorough  experimental  study  is  necessary  to 
validate  our  mathematical  model.  This  experimented  study  should  define  more  precisely  the 
chemical  nature  of  the  adhesion,  and  examine  a  wider  variety  of  materials. 
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Gm  (MPa) 


Figure  39.  Rupture  length  Xr  as  a  function  of  Gm 
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Gm  (MPo) 


Figure  40.  Damage  length  Xd  as  a  function  of  Gm 
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Figure  41.  Ineffective  length  iut  function  of  G„ 
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VI.  SUMMARY  AND  CONCLUSION 


The  purpose  of  this  reseach  was  to  improve  the  understanding  of  the  fiber-matrix  adhesion 
in  fibrous  composite  materials,  and  to  develop  a  micromechanical  model  of  the  matrix-fiber  stress 
transfer.  Attention  was  focused  on  the  stress  concentration  near  fiber  breaks  and  on  the  damage 
in  the  matrix  and  in  the  interphase. 

In  order  to  achieve  these  goals,  the  research  program  was  designed  by  combining:  (1)  two 
models  of  matrix-fiber  stress  transfer,  applicable  to  two  different  types  of  fibrous  composites,  (2) 
mechanical  tests  on  single  fiber  samples,  (3)  in  situ  optical  microscopy  under  polarized  light,  and 
(4)  elaboration  and  mechanical  testing  of  a  new  type  of  sample  containing  seven  fibers. 

In  the  analytical  part  of  this  research,  a  mathematical  model  was  developed  that  includes  the 
mechanical  properties  of  the  interphase,  the  stress  concentration  near  fiber  breaks,  and  the  elastic- 
plastic  behavior  of  both  the  matrix  and  the  interphase.  The  concept  of  ineffective  length  mentioned 
in  this  section  led  to  the  series  of  measurements  taken  under  polarized  light.  These  measurements 
appeared  to  correlate  better  with  the  model  than  did  the  critical  fiber  length.  A  second  model  was 
then  studied  to  take  into  account  the  slippage  of  the  end  of  the  fibers  in  the  matrix.  T  his  scheme, 
based  on  elasticity  theory,  is  only  valid  for  a  fiber  volume  fraction  small  enough  so  that  the  matrix 
shrinks  around  the  fibers.  This  is  the  case  for  a  single  fiber  specimen,  but  it  is  not  true  for  most 
of  the  fibrous  composite  materials  in  which  the  fibers  often  touch  each  other.  Such  a  model  must 
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be  elaborated  in  order  to  be  combined  with  our  "elasto-damaged"  model.  A  combination  of  these 
two  theories  would  constitute  a  fairly  complete  treatment  of  the  matrix-fiber  stress  transfer. 


The  results  of  the  critical  fiber  length  tests  that  were  performed  on  single  fiber  samples  dem¬ 
onstrate  that  there  is  no  relationship  between  this  test  and  the  matrix-fiber  adhesion.  The  results 
gave  only  quantitative  information  on  the  matrix-fiber  stress  transfer.  The  same  tests  run  on  the 
seven  fiber  specimens  were  even  less  significant  because  of  teh  premature  fracture  of  the  sample  due 
to  multiple  coplanar  fiber  breaks.  This  type  of  multiple  fiber  specimen,  however,  probably  has  a 
promising  future  because  of  the  potential  information  it  contains  about  fiber-fiber  interaction  and 
matrix  cracking. 

The  observation  under  polarized  light,  of  the  matrix  near  the  fiber  breaks,  gave  precise  and 
reproducible  results  that  correlate  with  the  predicted  value  given  by  the  "elasto-damaged"  model. 
This  observation,  performed  with  a  high  quality  optical  microscope,  is  the  most  accurate  means  of 
investigating  the  matrix-fiber  stress  transfer.  This  is  not,  however,  a  universal  test  method  because 
it  requires  that  the  matrix  be  a  birefringent  material,  which  is  not  the  case  in  most  composites. 
Consequently,  along  with  a  thorough  chemical  investigation  of  the  fiber-matrix  adhesion,  a  me¬ 
chanical  test  method  needs  to  be  developed.  This  new  test  should  involve  the  observation  of  a 
precisely  defined  phenomenon,  and  should  also  be  applicable  to  the  variety  of  materials  commonly 
found  in  commercial  composites. 
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APPENDIX  I. 


CHANGE  OF  COORDINA  TE  SYSTEM  FOR  THE 
MODELS  OF  THE  STRESS  FIELD  IN  THE 
MA  TRIXf  THE  INTERFACE ,  AND  THE  FIBER , 
DEVELOPED  BY  P.  J.  S 


Figures  A1  and  A2,  respectively,  represent  P.  J  Sabat's  coordinate  system  and  the  one  used 
in  this  report.  The  following  mathematical  manipulation  can  be  described  as  a  change  of  origin  and 
an  inversion  of  the  X  axis. 

Consequently,  the  equations,  p.  21  and  22  of  reference  6,  because  they  describe  the  state  of 
equilibrium,  are  not  affected  by  this  change  of  coordinate  system.  On  the  other  hand,  the  differ¬ 
ential  equation 


dl  rm 
dx1 


2 

-n  T« 


=  0 


(1/1) 
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where 


j2  2  RJEa  {R2a-R2m)  +  2  mfRf 
n  (Rm-Ri)lGm  +  p(Rl-Rf)IGl 


(2/1) 


has  to  be  solved  according  to  a  precise  coordinate  system  because  of  the  boundary  conditions.  The 
solution  of  Eq.  1 A  is  of  the  form 


Tm(*)  =  A  exp(  -  qx)  4-  B  exp(>/x) 


(3/1) 


and  the  boundary  conditions  are: 


-  Ef 


lim  rm(x)  =  0  ;  <7/(0)  =  0  ;  lim  o Ax)  =  o  — 7- 

X  — *  oo  *  J  r-*cm  J 


The  first  boundary  condition  and  eq.  3 A  imply: 


rm(x)  =  A  exp(  —  tjx) 


(4/1) 


The  equation  representing  the  equilibrium  of  the  fiber  is  now  recalled: 


f>  -r  ,  tG  do(X )  n 

2rr  Kr  1 1  +  n  Hr  — - - =  0 

J  J  ax 


(5  A) 


P.  J.  Sabat  also  assumes  that 


TiW  =  P  Tm(x) 


(6A) 


Including  eq.  6 A  into  eq  5 A  results  in: 


d  «;{x)  _  2/? 


dx  R,  m 


Tm(*) 


(7/1) 


111  sorting  eq.  4 A  into  eq.  7 A  and  integrating  gives: 


APPENDIX  I. 


104 


=  A(exp  [->/jc]  +  /0 


The  second  boundary  condition  along  with  eq.  8A  leads  to: 


«j(x)  =  ( exPC  ~  *1*]  ~  1) 


Inserting  the  third  boundary  condition  in  eq.  9A  gives: 


_  £/  -2 PA 

°  Ea~  Rfl 


—  a  Ej-Rj-y] 

2  VEa 


Finally  the  stress  in  the  fiber  and  in  the  surrounding  matrix  can  be  described  by 


a/x)  =  ( 1  -  exp[  -  */*]) 


—  a  RfErf\ 

TmW  - 2 ~JE -  CXPt-  ~  ^ 


where  rj  takes  the  value  expressed  in  eq.  2A. 
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